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We study the statistics of the relative separation between two fluid particles in a spatially smooth 
and temporally random flow. The Lagrangian strain is modelled by a telegraph noise, which is a 
stationary random Markov process that can only take two values with known transition probabilities. 
The simplicity of the model enables us to write closed equations for the inter-particle distance in 
the presence of a finite-correlated noise. In ID, we are able to find analytically the long-time growth 
rates of the distance moments and the senior Lyapunov exponent, which consistently turns out to 
be negative. We also find the exact expression for the Cramer function and show that it satisfies 
the fluctuation relation (for the probability of positive and negative entropy production) despite 
the time irreversibility of the strain statistics. For the 2D incompressible isotropic case, we obtain 
the Lyapunov exponent (positive) and the asymptotic growth rates of the moments in two opposite 
j^jT)' limits of fast and slow strain. The quasi-deterministic limit (of slow strain) turns out to be singular, 

i—{ , while a perfect agreement is found with the already-known delta-correlated case. 

■ PACS numbers: 47.52.+j, 47.27.-i, 47.27.cb 
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I. INTRODUCTION 



Understanding the statistics of the relative separation between two fluid particles in a generic random flow is a 
difficult task of significant importance [3, H, 0, 0, IE]- Here we treat small distances between particles, where the flow 
can be considered spatially smooth (so-called Batchelor regime @). We thus consider a general problem of dynamical 
systems where the linearized equation for the (infinitesimal) relative separation between two fluid particles, i?(t), 
takes the form 

> . 

■ d t R(t) = a(t)R(t) . (1) 

Here the matrix of the velocity derivatives o~(t) is called strain, it is taken in the co-moving (Lagrangian) frame and 
its statistics is supposed to be given. At times much exceeding the correlation time of the strain, general properties 
of systems with a dynamical chaos (called Lagrangian chaos in fluid mechanics) are described by the multi-fractal 

■ spectrum 7(n) of the growth rates of the moments: 

(R n (t)) ~ e^™)* . (2) 
Ch ■ The function 7(71) is convex, its derivative at zero is called the senior Lyapunov exponent, 

A = lim [(In R{t))/t] = 7 '(0) . (3) 



t — >oo 



Apart from that, not much is known about the general properties of the function 7(71). This is why few particular 
solvable cases are of much importance [3,0]- For the statistically isotropic strain, only three cases were found amenable 
to an analytic treatment [J]: the completely solvable short-correlated case in arbitrary dimension (where "f(n) is a 
parabola), and two partially solvable cases, the long-correlated case in two dimensions (where the Lyapunov exponent 
and its dispersion are found @) and the large-dimensionality case (where the second and the fourth moments are 
found together with the different-time correlation functions 0]). 

In this paper we describe another situation which allows exact analytical calculations: the telegraph-noise model 
for the strain [10(. This model enables one to describe analytically the dependence of j(n) on the correlation time 
of the strain. We shall study both the ID (compressible) and the incompressible 2D cases. The telegraph noise is 
a stationary random process, a(t), which satisfies the Markov property and only takes two values, a\ and 02- The 
probability, per unit time, of passing from the latter state to the former (or viceversa) is given by v\ (yi, respectively). 
In what follows, we shall consider some special cases, for which simple formulas hold. 
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(i) If ai = — a-2, then the process itself is stochastic but its square is deterministic, keeping the constant value 
a 2 (t) = a 2 (with a = \ai\ = | a.2 1 ) - We shall denote such processes with a star, e.g. a*(t). 

(ii) If the average (a) — (viai + ^2^2)/^ vanishes (with v = v\ + ^2), then the autocorrelation takes the form 
(a(t)a(t')) = e - "'* - ' '(^iaf + v^a^jv. For any analytical functional F[a] (function of time), a simple "formula of 
differentiation" [l(| then holds in this case: 

d t (a(t)F[a}) = (a(t)d t F[a\) - u{a(t)F[a]) . (4) 

We shall denote such processes with a tilde, e.g. a(t). Note that it is always possible to reduce to this case from a 
generic telegraph noise (in particular from one a*(t)), simply by subtracting from the latter its mean value. 

(iii) If both the properties in (i) and (ii) hold simultaneously, the maximum level of simplification is reached. For 
the sake of simplicity, we shall denote such processes not with a*(t) but rather with £(i). 

This type of noise has often been used to mimic more general colored noise in stochastic theories and has also 
been applied, e.g., to optics (see [TT[ | and references therein). An interesting application of the "telegraph concept" to 
turbulence can be found in [12| ■ In [10] a possible extension to the case of a multi- valued telegraph process is discussed. 
While the telegraph model is still far from any realistic flow, it is one step closer to reality than the previous (short- 
correlated) models that allowed for an analytical treatment of the two-point separation (4J. Our results provide, in 
our opinion, a first step in the analytic description of the role played by the strain correlation time in the statistics of 
the fluid-particle pair separation. 

The paper is organized as follows: in section [TT] we will study the ID case. In section IIIII we will move to 2D, 
describing the general incompressible isotropic case; the latter will be analyzed at first in its complete form, then 
(respectively in subsections MI AI and IIIIBI) by means of a simplified dynamics and eventually in the short-correlated 
limit. The quasi-deterministic limit is also shown in subsection IIII Cl Conclusions follow in section ITVl 

II. ID CASE 

The simplest meaningful model, making use of the telegraph noise, for the ID case, consists in assuming a strain 
function with the property (i), i.e. 

a{t) = a*(t). (5) 

Indeed, the problem is compressible locally, i.e. within every single realization of the noise, but for fluid-mechanical 
applications we want the average (ID) "volume" (R(t)) to be conserved. This means that we have to impose the 
constraint 

7 (1) =0. (6) 

From (fTJ) and ([5]) , it is straightforward to obtain the system of the two first-order ordinary differential equations which 
describe the time evolution of moments. It is indeed sufficient to multiply the starting equation by the noise itself 
(taken at the same time) in order to get a closed system, thanks to the simplicity of the telegraph form. Namely, if 
we define «o = (a*(t)} and rewrite accordingly ex.* if) = a(t) + olq, we obtain 

d t R n (t) = na*(t)R n (t) = n[a(t) + a Q ]R n {t) (7) 

and thus 

d t (R n (t)) = n(a(t)R n (t)) + na (R n (t)) . (8) 

We must then look for the evolution equation of (a(t) R n (t)) , which is readily done by exploiting ((4]): 

d t (&(t)R n (t)) = (a(t)d t R n (t)) - u{a{t)R n {t)) . (9) 

Taking ([7]) into account, we get (a(t)dtR n (t)) — na 2 (R n (t)) — no%{R n {t)) — nao(a(t)R n (t)), where we made use of 
the property stated in (i). Consequently, ^ is rewritten as 

d t (a(t)R n (t)) = n(a 2 - a 2 )(R n (t)) - (y + na ){a{t)R n {t)) . (10) 
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Equations ([8]) and (|T0|) constitute the system we were looking for. It is easy to recast the latter as a second-order 
differential equation for the quantity (R n (t)) alone, which can then be solved exactly: 

df (R n (t)) + ud t (R n (t)) - n(na 2 + va ){R n {t)) = , 
(R n (t)} = b n exp {-( ^Jv 2 + An 2 a 2 + Anva + v)t/2\ 

+/3„cxp|( v / ^ 2 + An 2 a 2 + Anva - v)t/2^ , 
-y(n) = (y/v 2 + 4n 2 a 2 + Anva — v)/2. (11) 
Constraint © applied to (|11[) implies 

a 2 ^Jv 2 + 4a 2 n(n ~ 1) - v 
a a = -— =S> 7 (n) = . (12) 

This is an even function of n — 1/2, it is convex and the asymptotics are 7(n)~ ± an at n — ► ±oo. 

The mean value «o is negative but its modulus cannot exceed a, which means that only situations with ajv (the 
analogue of the Kubo number) smaller than unity are physically relevant. Another interesting consequence is the 
relation v\ — V2 = —a, which reduces from three to two the number of free parameters in the original definition ([3]) 
of the strain. This means that the lower level of the noise persists longer: in the Lagrangian frame, indeed, particles 
spend longer time in contracting regions than in expanding ones, thus enhancing the weight of the former in the 
statistics. Accordingly, the Lyapunov exponent, computed as the derivative of (|12p taken in the origin, turns out to 
be negative as is required by the general theory [H, [3] : 

\ = i{Q) = -a 2 /v . (13) 

The joint probability density function P(R, t) can also be investigated. It satisfies the second-order partial differ- 
ential equation (Io| 

d 2 t P + vd t P - a 2 d R {Rd R {RP)] + a 2 d R (RP) = . (14) 

Equation (114j) has no stationary solutions, as both P = const, and P oc R^ 1 are non-normalizable. The asymptotic 
solution has a large-deviation form 

PfJI.il^e-^, (15) 

where X = t~ Y \n[R/ R(0)] can only belong (because of (|TJ) ) to the interval [—a, a]. The Cramer function H(X) is 
simply the Legendre transform of j(n): 



, > x_VEE*Wim. (16) 

2 2 la 

It is a convex function, which exactly satisfies the Evans-Searles-Gallavotti-Cohen fluctuation relation H(X) — 
H(—X) = X flU [IB] despite the fact that the strain statistics is not time- reversible [13] • The function vanishes 

quadratically at the minimum, represented by the Lyapunov exponent, H(X) X ^ X (X — \) 2 i/ 3 /4a 2 (is 2 — a 2 ), and it 
approaches vertically the boundaries of its compact domain: 



tti v \ x~±a v± a v 2 - a 2 i — 

Plugging (jTBJ) into (JT5J) , we obtain the final result: 

P(R, t) |i?|- 1 / 2 exp | (J{1 - a 2 jv 2 ){\ - \n 2 [R/ R{0)]/a 2 t 2 } - lj i^/2 J . 

It is worth mentioning what happens if one does not impose the constraint ((6]) and assumed a strain cr(t) = £(t) 
with the property described in (iii). In this case, 7(n) = (V v 2 + 4n 2 a 2 — v)/2 and the Lyapunov exponent would 
vanish. The Cramer function would be H(X) = v/2 — v\/ a 2 — X 2 /2a, so that 

P(R,t) liZj- 1 exp 



\[l -ln 2 [R/R(0)]/a 2 tA vt/2 



would be the solution to the equation for the joint probability distribution, similar to (|14p but without the last term. 
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III. 2D CASE 



In the 2D case, it is possible to consistently impose the incompressibility constraint locally, for every single realization 
of the noise. Since this case is most interesting for applications, we shall only consider this situation, which implies 
the strain matrix to be traceless. 

Consider a general statistically isotropic flow described by the strain matrix 



/(*) 



The noises £n(t) and (,m(t) are independent of one another but share the properties in (iii) with the same 

coefficients a and v, which can now range from to oo independently. Differently from the ID case, we will see that 
this is enough to mimic a realistic situation. The matrix a thus turns out to be singular and nilpotent: det a = = a 2 . 
By means of simple manipulations, analogue to those described for the ID case, one can write down a closed system 
for the moments of the components r n ^(t) = R\(t)R^ k (t) at every n, in the form 



d t (n 



{n) 



(17) 



The column vector VP^> has dimension 8(n + 1), being made up of: the n + 1 components of r n ,, itself, the 3(n + 1) 
components of r ni . times one noise, the 3(n + 1) components of r n ^ times two (distinct) noises, and finally the n + 1 
components of r„ i# times all the three noises. No other quantity is required to close the system, because as soon as 
a noise appears as square it becomes a deterministic quantity and can be taken out of the statistical average. The 
matrix A^") is thus of order 8(n + 1) x 8(n + 1) and turns out to be dependent only on a and v. The number of noises 
sharply increases in higher dimensions, thus strongly enhancing the number of components in 7ZS n ); this is why we 
confine ourselves to the 2D case. 

To study the long-time evolution of (r„,fc), one should identify the eigenvalue of with the largest positive real 
part. This is not an easy task from the numerical point of view for growing n, and it does not look feasible analytically 
for a generic n, even if the matrix itself can be easily written down as a function of n (the details can be found in [18| 
where, in particular, the whole spectrum of eigenvalues is found for n = 1 and n — 2). 

The evolution of the linear components (Ri) and (R2) is straightforward to describe: the largest eigenvalue is 0, 
which means that they do not show any exponential growth at large times; however, as they can change sign, this 
does not yield any information on (i?), i.e. on 7(1). The situation is more interesting for the quadratic components 
(R\ , -R1-R2, ^i); because in this case the largest eigenvalue is given by 



l(2) = -v+- 



{/2l6a 2 v + 3V5184aV 2_ 3^ 



y/21Qa 2 v + 3V5184a 4 ^ 2 ^Tl 



(18) 



which is always positive (i.e. describes an exponential growth of (R 2 ) = (R\) + (i? 2 .)) and behaves asymptotically as 



7(2) 



a" 
I — 

v 



96- 



2\/2a 2 v 





-0 




1 - 
















1 - 













2/3" 



for 
for 



(19) 



A. Exact reduced dynamics 

Since the general process is computationally very demanding, we succeeded to find a reduced set of variables 
invariant with respect to the evolution, i.e. a subset of components of IZ^ such that the corresponding rows in A^™) 
are nonzero only in columns whose index has been taken into account in the reduced subset. Consider even n: 

iT = (R 2 + RlT' 2 = £ (20) 



by definition for positive even n, thus the values 7(71.) can be extracted from the knowledge of n/2 + 1 of the first 
n + 1 components in TZ^ . On the contrary, this is not the case for the odd n's, which behave similarly to non-natural 
n's because they involve nonlinear operations (in this case a square root) that do not commute with the statistical 
averaging. 
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In particular, we analyzed the values n — 4, n — 6 and n = 8. Adopting a technique described in [Ta |. one 
gets matrices of order 5, 7 and 9 respectively (i.e. n + 1), whose eigenvalues with largest positive real part behave 
asymptotically as: 



7 (4) 



2 4 

a a a 

24 576^ for - -> 

^ 

a 



2V36a^ for 



7(6) 



48^! _ 2016^ for- 

V V 4 V 

a 
v 



2V1800a6i; for 



7(8) 



80— - 5280^ for- 

V V 4 V 

a 
v 



2v / 181440a i V for 



These relations provide us with three points of the curve j(n) in the asymptotic conditions, to which three more can 
be added: besides the already-mentioned 7(2), one should indeed remember that both 7(0) and 7(— 2) (in 2D) must 
vanish 0, 0, El- A simple extrapolation then suggests the following asymptotic behaviour: 

G(n)-+£(n)4 for-^0 
7(71) ~ { v v v (21) 

. g (nK l /(" +1 V/(™ +1 ) for £ -> 00 . 

Note that the noise correlation time is T = l/i>. The n-th growth rate grows linearly with T when the latter is small 
and decays as T _1 /' ,l+1 ^ when it is large. 

The coefficient G{n) can be proven rigorously to be n(n + 2) (see the next subsection), while such a proof does 
not look feasible for Q(n). However, starting from the quadratic form of the former, one can guess a fourth-order 
polynomial for the latter, Q{ri) = 3n(n + 2)(n 2 + 2n + 8)/4, and the correctness of this expression seems to be 
confirmed by the fact that it satisfies the six conditions in n = —2,0,2,4,6,8 despite possessing only five degrees of 
freedom. As the following coefficients in the power expansion at small ajv are expected to show higher powers of n, 
such an expansion is not uniform in the sense that it must fail for some (large enough) value of n. On the contrary, 
the coefficient g(n) must originate from some combinatorics, but has not been identified yet, while the presence of 
non-integer powers of a and v suggests that the limit a/v — > 00 is singular. From (|2ip . exploiting the knowledge of 
G(n) and the vanishing of g{0), we get the asymptotic behaviour of the Lyapunov exponent: 

2 a - for^O 

A = I v v (22) 



g\0)u for 



a 



B. Short-correlated case 



The results obtained in the previous subsection are exact, but they have been obtained for particular values of n, 
and the generalization of the coefficients (|2ip as functions of n was empirical. There is however at least one limit 
in which rigorous analytical results can be obtained: it is provided by the short-correlated case, corresponding to 
small ajv. We remind that in the exactly delta-correlated case 0, Q the probability distribution of R(t) is lognormal 
and the Cramer function is quadratic. Note that the correct limit is such that a/v —>■ but a 2 jv is finite. The 
key point in this limit is that the leading behaviour can be extracted rigorously, for any natural n, from a restricted 
dynamics, namely taking into account only the first 4(n+ 1) components in TZ^ (corresponding to the average of the 
coordinates r„ fc, and of the latter times one noise). The details of the calculations can be found in [18]; here we only 
mention that, upon introducing the small-Kubo Taylor expansion, we face a degenerate perturbation theory. Indeed, 
the eigenvalues corresponding to a fixed n are not definitely separate for a/v — > 0, but rather some of them vanish in 
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this limit, and we have to identify the one which keeps the largest positive real part in this process. The result of this 
analysis is that in the evolution equation for (r n one should account for the complete form, while for (£«r"n,fc) one 
can neglect all source terms involving two different noises: the system thus closes at this stage, reducing the order 
of the corresponding matrix from the initial 8(n + 1) to 4(n + 1). In other words, only the upper-left quarter of the 
matrix A^") is relevant. The problem can now easily be recast as a system of n+ 1 second-order differential equations 
for r n ^k- Indeed, we have 

d 2 r n , k = kR k 1 - 1 R^- k d 2 t R l + {n-k)R k 1 R^- k - 1 A 2 t R 2 

+k(k - l) J R^- 2 i?™- fe (d tJ R 1 ) 2 + (n - fc)(n — k- l)R k R^ k - 2 {d t R 2 ) 2 
+2k(n - k)R k - 1 R^- k - 1 d t R 1 d t R 2 , 
and, keeping into account the aforementioned simplification in the equations for (£«r rai fc), 

d 2 (r n ,k) + ^d t (r„ jfc ) = a 2 [(n 2 + 6k 2 - 6nk - n)(r n ^ k ) (23) 
+3fc(fc - l)(r„,fc_2) + 3(ra - k)(n - k - l)(r n ,k+2)] 
EE a 2 M^,(r r , k ,). 

The matrix M^ n \ defined by the right-hand side of (|2"3")l . is a priori of order (n + 1) x (n+ 1). However, it is apparent 
from (|23p that, for a fixed n, the dynamics of the components of (r n ,k) with even k are independent of those with odd 
k. Therefore, as our main goal is to reconstruct objects like (|20|) . we can just focus on subset M^™) of IVK") consisting 
of the lines with even k, and thus reduce to a (n/2 + 1) x (n/2 + 1) problem, simply neglecting the lines corresponding 
to odd k. In other words, we rewrite (|2"3"|) as 

d?(r„, 2fe } + vd t {r n , 2k ) =a 2 £ Mg fc ,(r n , fe ,) = a 2 M^(r„, 2fe ,) . (24) 

fe'=0,2,...,n 

Moreover, it is easy to show that the largest-positive-real-part eigenvalue \i n (divided by a 2 /v) that we are looking- 
for must belong to the spectrum of the subset matrix M W . Such subset matrix is tri-diagonal central-symmetric but 
not symmetric, which means that left and right eigenvectors differ. It is a simple task to prove that the value n(n + 2) 
(guessed for /j, 2 /(a 2 jv) asymptotically from the previous subsection) is actually an eigenvalue, but it is not possible 
to use a variational method to show that it has the largest positive real part. However, this is not needed, because 
the row vector built with the n/2 + 1 components of the binomial coefficient ( n [ 2 ) (for k = 0, 1, . . . , n/2) happens to 
be the left eigenvector of M^ n ) corresponding to the eigenvalue n(n + 2). Therefore, using (|20|) and (|24| . one has (for 
details see |18|): 

d 2 (R n ) + vd t (R") = n(n + 2)a 2 (R n ) . 

This implies 

as expected. Notice that the next-leading correction cannot be captured at this stage, because of the approximation 
introduced previously. To take it into account correctly, one should reformulate the analysis performed in this section, 
including also quantities like (^•^o^n,fc) but however excluding (£i£ii£uir nt k) 5 i-e. considering a reduced dynamics in 
the first 7(n + 1) components of TZ^ n \ 

It is worth mentioning that an analysis like the previous one, but applied to odd n's, leads to 

2 

a/i/->-0 . . . ,a 

[i n - (n— l)(n + 3) — , 
v 

in accordance e.g. with the exact result [i\ — found with the complete dynamics. However, as already pointed out, 
such values are not related to the curve "f(n). 

C. Deterministic limit in 2D 

The deterministic case corresponds to v = 0, when the noise £(i) takes a constant value ±a. In 2D, one gets 
d 2 R\ = = d 2 i? 2 (for any of the 2 3 = 8 possible combinations of the noise signs), thus both components are linear 



pin 



n(n + 2) — 
v 
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in time and R n t ^°° i". Consequently, a power-law temporal dependence is found for the separation and A = 0. It is 
worth noticing that this result is due to the fact that, with our choice, a 2 = 0, and is characteristic of the 2D isotropic 
situation. Therefore the limit is singular: an exponential growth is expected in general for finite a/v, but a power law 
is found when this ratio is infinite. As a result, perturbation theory and expansion in integer powers of v /a does not 
work. In (|2ip this is reflected by the presence of non-integer powers and the difference between the limits a — > oo and 



IV. CONCLUSIONS 

By assuming the telegraph- noise model for the velocity gradient (or strain matrix), we were able to carry out 
analytical computations and to obtain several results on the separation between two fluid particles. Focusing on 
smooth flows (Batchelor regime), wc firstly analysed the one-dimensional compressible case, finding explicit expressions 
for the long-time evolution of the interparticle-distance moments (|12p . the Lyapunov exponent (|13[) and the Cramer 
function ()16|) . Then we concentrated on the two-dimensional incompressible isotropic case, where a thorough analysis 
on the complete dynamics was made for the evolution of linear (n = 1) and quadratic (n = 2) components. Due to 
high computational cost, at higher n we focused on a restricted, though exact, dynamics: however, only specific values 
of n could be studied in this way, leading to the extrapolations (|21[) and (|22p . Such guess was rigorously proved in 
the quasi-delta-correlated limit, for which approximated equations were introduced. 

We believe that the present paper represents an interesting example of the use of a coloured noise for which the 
closure problem can be solved analytically. In particular, we established the dependence of the distance growth rates 
and of the Lyapunov exponents on the correlation time. 
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